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Сегментация примерно однородных 
по яркости областей цифровых изображений” 


Исследованы алгоритмы двунаправленной фильтрации, фильтр Кувахары, алгоритм среднего сдвига для 
сегментации примерно однородных по яркости областей цифровых изображений. Представлены 
результаты экспериментов. Даны рекомендации по выбору параметров. Результаты сегментации могут 
использоваться в качестве входных данных для высокоуровневых процедур автоматической обработки. 


Введение 


Целью данной статьи является поиск и разработка методов и алгоритмов, позво- 
ляющих автоматически сегментировать изображения на примерно однородные по яр- 
кости области. Эти области помогают определить на изображении границы некоторых 
объектов. Такая сегментация является важным предварительным этапом при решении 
многих прикладных задач дешифрирования и распознавания изображений. 

Целью сегментации является разделение изображения 5 на визуально однородные 
области 5;. Степень однородности зависит от решаемой задачи. 


© = {5 },1=1, 
Ад] 9, 
ТЫ М 


УрЛЕЬМ:1 = Р(5, 15, ) = ложь. 


Сегментация большинства реальных изображений представляет собой сложную 
задачу. Многие существующие алгоритмы сегментации изображений основываются на 
одном из двух базовых свойств сигнала яркости: разрывности и однородности [1]. 
В первом случае подход состоит в разбиении изображения на основании резких изме- 
нений сигнала, таких как перепад яркости на изображении. Во втором — используют раз- 
биение изображения на области, однородные в смысле заранее выбранных критериев. 

В наших экспериментах в дополнение к широкоиспользуемым методам сегмента- 
ции исследовались алгоритмы двунаправленной (Бабега!) фильтрации [2], фильтр 
Кувахары [3], алгоритм среднего сдвига (теап $1) [4]. 


* Работа частично выполнялась в рамках проекта ПМТА$ № 06-1000024-9100. 
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1. Двунаправленная нелинейная фильтрация 


Двунаправленная фильтрация выполняет сглаживание объектов на изображении, 
при этом сохраняет структуру (края объектов) изображения. 

Двунаправленная фильтрация — это нелинейный и неитерационный процесс, 
комбинирующий пространственную (дотла) и яркостную (тапое) фильтрацию. Двунап- 
равленная фильтрация изображения /(х) в общем виде описывается выражением: 


оз, | 


где 1(х) — выходное изображение, с(&,х) и 5(1(2), Г(х)) — геометрическая близость 
и фотометрическое сходство (рВоютейлс зппЙагиу) соответственно между центром 
окрестности х и соседней точкой &. 

Двунаправленная фильтрация заменяет исходное значение каждого пикселя изоб- 
ражения на сглаженное значение, вычисленное в окрестности заданного размера. 

В простом случае двунаправленной фильтрации в качестве функции близости 
с(е,х) используется гауссиан: 


й(х 


с(о,х) =е | 
где 4 (Е ь х) =а (Е _ х) = |Е _ х| — евклидово расстояние между пикселями & их. 


Аналогично, в качестве функции сходства 5( Л (Е ), У (х)) используется гауссиан: 


(но 
$(о,х)=е о - 


где 6(9, /)=5(- Г)=|#- /| — евклидово расстояние между значениями яркостей ф и/ 


Алгоритм двунаправленной нелинейной фильтрации 


1. Значения яркости исходного полутонового изображения А размером МхМ№ 
привести к диапазону [0,1]. 

2. Задать размер маски фильтра у, допустимые среднеквадратичное отклонение 
оди среднеквадратичное отклонение яркости пикселей о,,. 

3. С помощью уравнения (1) сформировать маску размером (2%+1)*(2и-+1) для 
пространственной фильтрации, хе (-\%,\), уЕ (-и, м»). 


_ 2 +92 


б(х,у)=е №. (0 


4.В цикле для каждого пикселя изображения А(1,)), где 1<1<М,1< 7 <М, вы- 


ПолнНиТЬ: 
а) с помощью уравнения (2) сформировать матрицу Н отклонений яркости в 
окрестности Г размером (2%+1)*(2+1) текущего пикселя А(1,/): 

(бъу-Ае ру 

202 


Н(х,у)=е | , (2) 
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6) вычислить матрицу Е путем поэлементного умножения Н на С: 
Е(х, у) = Н(х,у)-С(х,у), (3) 


в) сформировать результирующее изображение В, вычислив в окрестности / 
размером (2*+1)*(2%-+1) каждого пикселя выражение (4): 


В(х,у)= 
(х, у) УР, 


5. Конец цикла. 
6. Конец алгоритма. 


В 


Рисунок 1 — Результат двунаправленной фильтрации изображения спутника Аз{ет: 
а, г- исходные изображения (компоненты 1 и 3 соответственно); 
6, д— изображения, обработанные фильтром, с параметрами: о, =2, а, = 0,01; 
в, е- изображения, обработанные фильтром с параметрами: о, =2, а, = 0,03 


На рис. | представлено исходное изображение спутника Азег и результаты 
двунаправленной фильтрации в окрестности 3х3 с различными значениями парамет- 
ров бу иИо,. Для демонстрации свойств фильтрации на исходное изображение нане- 


сены черные и белые точки. На рис. 16 - в видно, что резкие края даже у малых 
объектов не размываются, но небольшие яркостные отличия сглаживаются. Такая 
фильтрация полезна для последующей классификации объектов изображения, напри- 
мер, при интерактивном дешифрировании. 


334 «Искусственный интеллект» 32008 


Сегментация примерно однородных по яркости областей... 4М 


2. Фильтр Кувахары 


Этот оператор выполняет нелинейную фильтрацию изображений с сохране- 
нием резких краев. Результат похож на грубую сегментацию изображения. 


Алгоритм фильтрации Кувахары 


1. В цикле по всем пикселям исходного изображения, в окрестности каждого 
пикселя применить маску размером (2+ 1)*(2%+ 1). 

2. Окрестность разделить на 4 области, указанные цифрами 1, П, Ш, ГУ, 
размером ии’. 


Рисунок 2 — Окрестность пикселя Х 


3. Вычислить среднюю яркость и дисперсию в каждой из областей 1, П, Ш, ТУ. 

4. Присвоить центральному пикселю Х среднее значение яркости той области, 
дисперсия которой наименьшая. 

5. Конец алгоритма. 


Примеры обработки изображений фильтром Кувахары приведены на рис. 3—5. 
На обработанных изображениях хорошо видна сегментация на примерно однородные 
по яркости области с сохранением четких границ. Такая обработка полезна для при- 
мерной классификации полей и других земель, а также в качестве исходных данных 
для более точной интерактивной сегментации отдельных объектов. 


а 


Рисунок 3 — Результат фильтрации Кувахары 
а исходное изображение спутника Азег (компонент 3); б — изображение, 
обработанное фильтром с параметром и = 11 
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Рисунок 4 — Результат фильтрации Кувахары 
а-— исходное изображение спутника Азег (компонент 4); б — изображение, 
обработанное фильтром с параметром = 5 


Рисунок 5 — Результат фильтрации Кувахары вышеприведенных компонентов 3 и 4 
спутника Азег, обработанных три раза одним и тем же фильтром с параметром и = 3 


3. Сегментация методом среднего сдвига 


Данный адаптированный алгоритм описан в статьях [4], [5]. Основные идеи 
алгоритма описаны ниже. 

Предполагаем, что функция плотности вероятности р(х) р-мерного вектора 
признаков х — унимодальна. Это предположение используется для упрощения, в 
дальнейшем оно может быть опущено. Сфера 5х радиуса г, с центром в х содержит 
вектора признаков у таким образом, что || у-х||< г. Тогда математическое ожидание 
вектора 2 = у—х при заданных х и 5х: 

р) (5) 


и=Е[2| 5%] = Ооо О 
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Если сфера 5х достаточно мала, можно записать: 


РО? 5х) = р(х)Г,, (6) 
где Г;, =с-г” — объем этой сферы. Аппроксимация р(у) первого порядка 
2) = р(х)+ (7-х) УрО), (7) 


где Ур(х) - градиент функции плотности вероятности х. Тогда 


(у-х)(ф-х)” Ур(®) 
= ау. 8 
` | Г. 2) - 
Значение интеграла: 
__г_ Ур (9) 
р+2 р(х) 
ИЛИ 
вые ИИС (10) 
р+2 р(х) 


Алгоритм среднего сдвига 


1. Выбрать радиус г окна поиска. 

2. Выбрать начальное местоположение окна. 

3. Вычислить средний вектор сдвига, передвинуть окно поиска на этот вектор. 
4. Повторить до сходимости. 


Рисунок 6 — Результат сегментации методом среднего сдвига 
а— исходное изображение спутника ОшсКЫта; 6 — результат сегментации 
методом среднего сдвига 
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На рис. 6-7 представлены результаты сегментации методом среднего сдвига. 
Данный метод не применим для площадных объектов, обладающих явно выражен- 
ными текстурными свойствами, так как текстурные характеристики будут заменены 
усредненным значением яркости. В случае разных текстурных свойств, но одина- 
ковых средних значений яркости в окне анализа, эти объекты могут быть объеди- 
нены в один. 


В 


Рисунок 7 — Результат сегментации методом среднего сдвига 
а — исходное изображение спутника А$ег (компонент 4); 6 — г - результат 
сегментации методом среднего сдвига с разными параметрами 


Выводы 


Из-за наличия шумов различной природы на объектах изображения наблюда- 
ются локальные яркостные искажения, связанные с неоднородностью отражения. 
Описанные выше алгоритмы сохраняют структуру и форму объекта, обеспечивая визуа- 
льное улучшение разделения площадных объектов, отличающихся между собой по яр- 
кости. Алгоритмы предназначены для предобработки снимков, содержащих площадные 
объекты, примерно однородные по яркости. Применение таких алгоритмов полезно для 
интерактивных ГИС-технологий, использующих дешифрирование различных объектов. 
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